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Abstract. We review and generalize the recent progress in a soliton cellular 
automaton known as the periodic box-ball system. It has the extended affine 
Weyl group symmetry and admits the commuting transfer matrix method and 
the Bethe ansatz at q = 0. Explicit formulas are proposed for the dynamical 
period and the number of states characterized by conserved quantities. 



1. Introduction 

In ^3] j a class of periodic soliton cellular automata is introduced associated with 
non-exceptional quantum affine algebras. The dynamical period and a state counting- 
formula are proposed by the Bethe ansatz at q = |H|. I n this paper we review and 
generalize the results on the An case, where the associated automaton is known as 
the periodic box-ball system (141 The box-ball system was originally introduced 
on the infinite lattice without boundary |18l 117) . Here is a collision of two solitons 
with amplitudes 3 and 1 interchanging internal degrees of freedom with a phase shift: 

1114221111131111111111111111 ■ 
1111114221113111111111111111 ■ 
1111111114221311111111111111 ■ 
1111111111114232111111111111 ■ 
1111111111111121432111111111 
1111111111111112111432111111 ■ 
1111111111111111211111432111 

The system was identified 0] Ej with a solvable lattice model at q = 0, which 
led to a direct formulation by the crystal base theory |S] and generalizations to the 
soliton cellular automata with quantum group symmetry [7||H]. Here we develop the 
approach to the periodic case in jT^j further by combining the commuting transfer 
matrix method and the Bethe ansatz j2| at q = 0. 

In section |2 we formulate the most general periodic automaton for g n — in 
terms of the crystal theory A commuting family of time evolutions {Tj} is intro- 
duced as commuting transfer matrices under periodic boundary condition. The asso- 
ciated conserved quantities form an n-tuple of Young diagrams m = (to*- 1 ), . . . , mS n '), 
which we call the soliton content. 

In section|3|we invoke the Bethe ansatz at q = JT] to study the Bethe eigenvalue 



A[ r) relevant to Tf. The Bethe equat ion is linearized into the string center equation 

and the A| r ' is shown to be a root of unity. We also recall an explicit weight mul- 
tiplicity formula obtained by counting the off-diagonal solutions to the string center 
equation jll) . It is a version of the fermionic formula called the combinatorial com- 
pleteness of the string hypothesis at q = 0. These results are parameterized with the 
number of strings, which we call the string content. 

l 
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In section 0] two applications of the results in section are presented under the 
identification of the soliton and the string contents. First we relate the root of unity 
in the Bethe eigenvalue to the dynamical period of the periodic automaton 

(r) 

under the time evolution T ; . Second we connect each summand in the weight 
multiplicity formula |11| to the number of states characterized by conserved quantities. 

In JSj; similar results have been announced concerning the highest states. Our 
approach here is based on conserved quantities and covers a wider class of states 
without recourse to the combinatorial Bethe ansatz at q = 1 JUj . We expect parallel 
results in general g n . In fact all the essential claims in this paper make sense also for 
0„ = Dn and Eq? 8 . Our formulas fTTJ) and (fT3")l include the results in ^3 proved 

by a different approach as the case g n — with B = and I = oo. For the 

standard notation and facts in the crystal theory, we refer to El Ej • 

2. Periodic A)P automaton 

Let i? aj (1 < a< n,j G Z>i) be the crystal |5] of the Kirillov-Reshctikhin mod- 
ule wj a) over U q {A ( n ] ). Elements of B aj are labeled with semistandard tableaux 
on an a x j rectangular Young diagram with letters {l,2,...,n + l}. For exam- 
ple when n = 2, one has B 1 ' 1 = {1, 2, 3}, B 1 ' 2 = {11,12,13,22,23,33}, B 2 ' 2 = 

{ 2 2.2 3.3 3.2 3.3 3.3 3 } as sets. AS(B^)={C d b\bGB^,dG 
Z} denotes the afhne crystal. The combinatorial R is the isomorphism of affine crystals 
AR(B a ' j ) ® AS{B b ' k ) ^ AS(B b ' k ) <g) AS(B a ' j ) [H] . It has the form R(( d b®Cc) = 
( e+H c (g) ( d ~ H b, where H = H(b ® c) is the energy function. We normalize it so 
as to attain the maximum at H(u aJ <g> u b > k ) = 0, where u a > j G B QJ denotes the 
classically highest element. We set B = B ri ' 11 ® B r2 ' 12 ® • ■ ■ ® b tl ' 1l and write 
AS(B ri h ) (g) • • • (g) AS(B rL ' lL ) simply as AS(B). An element of B is called a state. 
Given a state p = &i <g> • • ■ ®&l G B, regard it as the element £°oi <g> • • • <g> £°ol G Aff(B) 
and seek an element v G B r < 1 such that ®P— (C dl b'i €5 • ■ ■ <8> ( dL b' L ) <g> C e ^ un- 
der the isomorphism Aff(B r ' i ) ® Aff(S) ~ Aff(B) <g> Aff(B r >'). If such a u exists and 
C, dl b' l <g) • • ■<S)C t d,L b' L is unique even if v is not unique, we say that p is (r, l)-evolvable and 

write 2j (p) = 6^ <g> • • • ® 6' L G -B and 25j (p) = e = — di di,. Otherwise we say 

that p is not (r, Z)-evolvable or T ; (r) (p) = 0. We formally set r/ r) (0) = 0. E\ r) G Z> 
holds under this normalization. Our automaton is a dynamical system on BU{0} 
equipped with the family of time evolutions {Tj | 1 < r < n, I £ Z>i}. is the 
q = analogue of the transfer matrix in solvable vertex models. It is invertible and 
weight preserving on B. Using the Yang-Baxter equation of the combinatorial i?, one 
can show (cf. 

Theorem 2.1. The commutativity T- 'T^ (p) — T^Tj a \p) is valid for any (a, j), (6, k) 
andp G B, where the both sides are either inBorO. In the former case Ej {T^ (p)) — 
E^{p) and E { k b) (T^ a \p)) = E { k b \p) hold. 

Thus, {E^ I 1 < a < n, j G Z>i} is a family of conserved quantities. 

Conjecture 2.1. For any l<a<n andp G B, there existsi > 1 such thatTj. (p) 7^ 
«/ and onZy i/ fc > i. T7ie Ziraii limfc^oo T k a \p) £ B exists and E\ a \p) < E^^p) < 
■ ■ ■ < Ej ^ (p) = Ej+i (p) = ■ ■ ■ holds for some j > i. 
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Let So, Si, . . . , S n be the Weyl group operators 5 and pr be the promotion operator 
acting on B component- wise. For instance for A^\ pr^ 334® 1344^ = 

4 4 4 ® 1124 € B 2 ' 3 (gi B 1 ' 4 . They act on B as the extended affine Weyl group 
W( Qn ) = W(A n 1) ) = (pi, So, Si,..., S n ). 

Theorem 2.2. IfT^ a) (p) ^ 0, then for any w E W(A n 1] ), the relations wT^ a) '(p) = 
Tj a \w(p)) and Ej(w(p)) = Ej(p) are valid. 

A state p G -B is called evolvable if it is (a, j)-evolvable for any (a,j). In Conjec- 
ture IOI we expect that the convergent limit equals a translation in W(An^). 
Compared with T; in |13| . the family {T^ } here is more general and enjoys a larger 
symmetry W(An^). Define the subset of B by 

(1) P(m) ~ {p E B \ p : evolvable, E^\p) = min(j, k)m^}. 

k>l 

Pictorially, m — (m^\ . . . ,m^) is the n-tuple of Young diagrams and E^ (resp. 
Tnfjf^) is the number of nodes in the first k columns (resp. number of length k rows) 
in nrS a \ We call m the soliton content. 

Remark 2.1. P(m) is W(A n ) -invariant due to Theorem \2.2L 

Given p E P(m), T^ a \p) E P(m) is not necessarily valid. For instance p = 
112233 e P(((22), (2))) C B = (B 1 - 1 )® 6 but t[ 2) (p) = 213213 is not evolvable since 
(T[ 2) ) 2 (p) = 0. On the other hand one can show 

Proposition 2.1. If p E P(m) and (rj o) )*(p) ^ for any t, then (T^ a) Y(p) E P(m) 
for any t. 

Let W, A a , a a be the Weyl group, the fundamental weights and the simple roots of 
A n , respectively. We specify p^ = Pj(m) by © and set 

n 

(2) A(m)=^pWA a , H = {(a,j)\l<a<n,jeZ>i,mf>0}. 

a=l 

Conjecture 2.2. P(m) 7^ if and only if p^ > for all (a,j) E H. {wtp \ p E 
P(m)} = WX(m). 

The claim on the weights is consistent with the VF^l^-invariance of P(m). 

Here is an example of time evolutions in B = B 1 ' 1 &B 1 ' 1 ^B^-^^B 1 ' 1 &B 1 ' 1 &B 1 ' 1 ® 
B 1 ' 2 with W(AP) symmetry. The leftmost column is p ,T^ ] (po),T {2) T^ ] (p ) and 
j,(S) j,(2) rpW^p^ from the top to the bottom. At each time step, the states connected 
by the Weyl group actions So and Si are shown, forming commutative diagrams. (• 
signifies ®.) All these states belong to P(((3111), (21), (1))). 

2 • 1 • 233 • 4 • 1 • 2 • 12 & 2 • 4 • 233 • 4 • 1 • 2 • 24 & 1 • 4 • 133 • 4 • 1 • 2 • 14 

1- 2- 123- 3-4- 1-22 4 • 2 • 234 • 3 • 4 • 1 • 22 4 • 1 • 134 • 3 • 4 • 1 • 12 
1 • 2 • 112 • 3 • 2 • 3 • 24 1 • 2 • 244 • 3 • 2 • 3 • 24 1 • 2 • 144 • 3 • 1 • 3 • 14 

2- 3- 112 -4 -2-3- 12 2 • 3 • 244 ■ 4 • 2 • 3 ■ 12 2 • 3 • 144 • 4 • 1 • 3 • 11 
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Remark 2.2. Let Rj(l <j< L — l) be the combinatorial R that exchanges the j-th 
and (J+l)-st components in B = B ri ' 11 ® • • • eg) B Tl ' 1l and 7r(6i ® • • • ® 6x) = 6x (g) 61 Cg> 
• • • ® 6x1-1. Together with Rq = ir~ 1 Riir, they act on U se e i .B rs i''< , i ® ••• ® B Ts l ,1s l 
as the extended affine Weyl group W(A^_ l ) = (71", Rq, . . . , Rl-i)- Theorem \2.°A is 
actually valid for w G W(g n =An) x ^(A^ 1 ^) as in [jj]. In i/ie homogeneous case 

(?*i,Zi) = ••• = (rx,^x), fie W^-A^ 1 ^) symmetry shrinks down to the ■K-symmetry, 
which is the origin of the adjective "periodic". 

3. Bethe ansatz at q = 

Eigenvalues of row transfer matrices in trigonometric vertex models are given by 
the analytic Bethe ansatz ^JE]- Qr{ u ) — rifc smn7r ( M — V~^ u k) De Baxter's 



"k 

Q-function where {v,j} satisfy the Bethe equation eq.(2.1) in [TT]. We set q = e~ 2lTh 
and C = e 27rtl . For the string solution ( 11 Definition 2.3), the relevant quantity to 
our r/ r) is the top term of the eigenvalue A ( f \u) (cf. 12 (2.12)): 

jo 

Here 2^ is the center of the a-th string having color a and length j. Denote by rrij 

the number of such strings. We call the data m = (rtij ) the string content. The 

product in @ is taken over j £ Z>i and 1 < a < rrij. -Ej is given by the same 
expression as in (JIJ as the function of m. At g = the Bethe equation becomes the 
string center equation (JT] (2.36)): 

(4) n n^S)- 4 — =(-i)^ ,+m " )+i , 

(6,fe)eif 0=1 

(5) A ajoit bk0 = SabSjkSapiPj + m^) + C ab min(j, fe) - 5 a b5jk, 

L 

(6) = ^min(j,li)6 ari - ^ C q6 mm(j, k)m k b) , 

i=l (b,k)£H 

where (C Q t,)i< a ,(,<„ is the Cartan matrix of j4„. To avoid a notational complexity we 
temporally abbreviate the triple indices aja to j, bk(3 to fc and accordingly to z/. 
etc. Then Q and @ read 

(7) a« = n(-**r> n(-^)^ fe = (- 1 ) 53 ' 

where pk is given by pfc = Sb r min(fc, I) for k corresponding to bk[3, and Sj is an integer. 
Note that Aj,k — Akj- Suppose that the q = eigenvalue satisfies (Aj*^)' p i ' = ±1 for 
generic solutions to the string center equation . It means that there exist integers £j 
such that J2j £jAj,k = V\ pk, or equivalently £j = d ^ c fj^ , where A[j] denotes the 
matrix A — (Aj,k) with its j-th column replaced by (pi, P2, ■ ■ ■ )■ In view of the con- 
dition V*y € Z, the minimum integer allowed for is v\ r) = LCM^l, \J k ' d ^ fc] 
where LCM stands for the least common multiple and U' k means the union over those 



V\ ' here should not be confused with the symbol in IrTl , 
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k such that A[k] 7^ 0. Back in the original indices, the determinants here can be 
simplified (cf. ^I] (3-9)) to those of matrices indexed with H: 



(8) V\ r) = UM\l, |J' 

(b,k)eB 



detF[6, k] 



(10) F[6,fc] 6 



where the matrix F = (F a j,bk)(aj),Q>> k )€ H is defined by 

(9) F aj , bfc = 5 a b5jkp\ a) + C ab min(j, . 

The matrix F[b, k] is obtained from F by replacing its (6, fe)-th column as 

_ \Faj, cm (c, m) 7^ (b, k), 

y S ar min(j, I) (c, m) = (6, k) . 

The union in (jSJ is taken over those (&, fc) such that det F[b, k] 7^ 0. 

The LCM © can further be simplified for A ( ^\r = 1, Vr t = VZ* = 1. We write 

just as and parameterize the set H = {j <E Z>i | > 0} as = {(0 <) J\ < 
■ ■ ■ < J s }. Setting if. — min( Jfc, I) and io = 0, one has 

(11) 7>, (1) = LCM I 1, M'^V , 

V ^0 

where < t < s — lis the maximum integer such that it+i > it. 

Let us turn to another Bethe ansatz result, the character formula called combina- 
torial completeness of the string hypothesis at q = [TTj : 

L 

(12) JJchB ri,ii = ^ft(m)e A(m) , 

2= 1 m 

1 (a) 1 v 

where (J) = s(s - 1) • • • (s - t + l)/t! and cb_B r >' is the character of B r - 1 . A(m), pj o) 

and F are defined by J3J, © and (|5J). The sum in 112(1 extends over all S Z>o 
canceling out exactly leaving the character of B. 112|) and (|13|) are the special cases 
of eq.(5.13) and eq.(4.1) in |TT], respectively. O(m) (denoted by R{v,N) therein) is 
the number of off-diagonal solutions to the string center equation with string content 



m. It is known (JJ Lemma 3.7) that f2(m) 6 Z>i provided that pj a ' ) > for all 
(a,j)&H. 



4. Dynamical period and state counting 

In fllfl. the soliton content m = (mj a ' ) ) is introduced as the conserved quantity 
associated with the commuting transfer matrices. One the other hand, the in 

the string content m = (m^ ) is the number of strings of color a and length j in the 
Bethe ansatz in section [21 From now on we identify them motivated by the factor 
C~ E ' in 121 and some investigation of Bethe vectors at q = 0. In view of Conjecture 
12.21 the data of the form m — (m^ ) is defined to be a content if and only if py' > 
for all (a,j) E H. Thus det F > and A(m) in © is a dominant weight for any 
content m. 
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Conjecture 4.1. If p € P(m) and (r/ r) )*(p) ^ for any t, the dynamical period of 
p under (minimum positive integer t such that (7^)*(p) = p) is equal to 



generically and its divisor otherwise. 



In the situation under consideration, the whole T, orbit of p belongs to P{m) 
due to Proposition 12.11 Naturally we expect (A^) 73 ! ' = 1, which can indeed be 

verified for Conjecture 14.11 has been checked, for example in case, for 

B = (i? 1 ' 1 )® 3 ® B 2 ' 2 and B 2 ' 1 <g> i? 2,1 ® B 3 ' 1 ® 5 3 > 2 . 

Let us present more evidence of Coniecture l4.ll To save the space, • = ® is dropped 
when B = (B 1,1 )® L . In each table, the period under Tj with maximum / is equal 
to that under !«, . 

1221121122221, content = ((321)) 



A^ 1 , state 



(r,0 


LCM of 


= period 


(1,1) 
(1,2) 
(1,3) 


1, 13, 13, 13 

i 91 91 91 
> 3 ' 16' 16 
i ni 273 273 


13 
91 
273 



A[ 1} , state = 122 • 112 • 12 ■ 1222 • 2 ■ 11111 • 1122 ■ 111, content = ((4321)) 



(r,0 


LCM of 


= period 


(1,1) 


1, 


2 






2 


(1,2) 


1, 


7, 


2 21 


42 


42 


(1,3) 


1, 


14, 


7 21 


21 
2 


42 


(1,4) 


1, 


21, 


21 63 
2 ' 8 ' 


126 
29 


126 



,(1) 



state = 134 • 34 • 1 



A 



(r,l) 




(1,1) 


1, 


(1,2) 


1, 


(2,1) 


1, 


(2,2) 


1, 


(2,3) 


1, 


(3,1) 


1, 


3 , state 


(r,0 




(1,1) 


1, 


(1,2) 


1, 


(1,3) 


1, 


(2,1) 


1, 


(2,2) 


1, 


(2,3) 


1, 


(3,2) 


1, 



134 • 23 • 1 • 13, content = ((432), (31), (1)) 
LCM of = period 



380 


95 


95 


380 


380 


380 


39 ' 


6 " 


6 > 


31 ' 


27 ' 


29 


190 


95 


95 


190 


190 


190 


39 ' 


12 ' 


12 ' 


31 ' 


27 ' 


29 


190 


95 


95 


190 


190 


190 


13 ' 


4 • 


4 ' 


137' 


9 ' 


73 


70 


38 


38 


76 


76 


70 


5 ' 


3 ' 


3 ' 


41 ' 


21 ' 


31 


95 


95 


95 


95 


95 


95 


6 ' 


11 ' 


11 ' 


34' 


48 ' 


41 


380 


95 


95 


380 


380 


380 


13 ' 


2 ' 


2 ' 


137 ' 


9 ' 


263 



380 
190 
190 
76 
95 
380 



= 233 • 



1 1 

3 4 



LCM of 


n 

2 ' 


11, 


22 


11 

4 ' 


11 

2 ' 


11 


11 


11 


22 


6 ' 


3 ' 


3 


11, 


33 


06 


7 ' 


7 


11 


33 


33 


2 ' 


14 ' 


7 


11 


11 


22 


3 ' 


7 ' 


7 


11, 


33 
7 ' 


33 
20 



• 1, content = ((3), (3), (2)) 
period 



22 
11 

22 
66 
33 
22 
33 
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Here, content =((3 111), (44), (2)) for example means that 
1, ml 2 '' = 2 and the other mi- s are 0. 



Q (1) 

= 3, = m 



,(3) 



Let us turn to another application of the Bethe ansatz results (|12|) and i|13|) . We 
introduce T(P(m)) = |J" =1 {J 3 >A T j a) (?) I P £ P(™)}, which is the subset of B 
consisting of all kinds of one step time evolutions of P(m). Under Conjecture 12.11 
any state p € P(m) is (a, j)-evolvable for j sufficiently large. Thus from Proposition 
12.11 p is expressed as p = {T^) k (p) for some k, showing that T{P{m)) 2 P(m). In 
general T(P(m)) can contain non-evolvable states which do not belong to P(m). 

Conjecture 4.2. For any content m such that T(P(m)) — P(m), the following 
relation holds: 



(14) 



n(m) = 



\P(m)\ 
\WA(m)\' 



In view of Remark |2. II and Conjecture 12.21 the right hand side is the number of 
states in the periodic automaton having the content m and a fixed weight. Thus 
it is equal to §{p S P(m) | wtp = A(m)}. In case T(P(m)) 2 P(m), we expect that 
|P(m)|/|FyA(m)| is a divisor of f2(m). 

Let us present two examples of Conjecture 14.21 In the periodic A3 automaton 
with B = B 1 ' 2 ® B 1 ' 1 ® B ia ® B 1 * 1 , there are 1600 states among which 824 are 
evolvable. They are classified according to the contents m in the following table. 



m 


A(m) 


\WX{m)\ 


\P(m)\ 


fl(m) 


(0,0,0) 


(6,0,0,0) 


4 


4 


1 


(W,0,0) 


(5,1,0,0) 


12 


48 


4 


((H), 0,0) 


(4,2,0,0) 


12 


24 


2 


((2), 0,0) 


(4,2,0,0) 


12 


72 


6 


((21), 0,0) 


(3,3,0,0) 


6 


24 


4 


((3), 0,0) 


(3,3,0,0) 


6 


36 


6 


((H),(i),0) 


(4,1,1,0) 


12 


96 


8 


((22), (2), 0)* 


(2,2,2,0) 


4 


24 


12 


((21),(1),0) 


(3,2,1,0) 


24 


432 


18 


((111), (11), (1)) 


(3,1,1,1) 


4 


16 


4 


((211), (11), (1)) 


(2,2,1,1) 


6 


48 


8 



In the second column, (Ai, A2, A3, A4) means A(m) = (Ai — A2)Ai + (A2 — As)A2 + 
(A 3 — A4)A 3 . In the last two cases, the subsets of P(m) having the dominant weight 
A(m) are given by 

{11-2-13-4, 12-3-14-1, 13-4-11-2, 14-M2-3} for m = ((111), (11), (1)), 
{11-2-23-4, 12-2-13-4, 12-3-244, 13-4-12-2, 
14-1-22-3, 22-3-14-1, 23-4-11-2, 24-1-12-3} for m = ((211), (11), (1)). 



In the case of B 
evolvable. 



B 2 ' 1 Cg) B 2 ' 1 ® B 2 ' 2 , there are 720 states among which 518 are 
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777, 


A(m) 


\W\(m)\ 


\P(m)\ 


n(m) 


(0,0,0) 


(4,4,0,0) 


6 


6 


1 


(0,(1), 0) 


(4,3,1,0) 


24 


72 


3 


(0,(n),(i)) 


(4,2,1,1) 


12 


36 


3 


(0,(2),0) 


(4,2,2,0) 


12 


48 


4 


((i),(n),0) 


(3,3,2,0) 


12 


36 


3 


((l), (11), (l)) 


(3,3,1,1) 


6 


72 


12 


(W,(21),(l)) 


(3,2,2,1) 


12 


240 


20 


((2), (22), (2))* 


(2,2,2,2) 


1 


8 


32 



The assumption T(P(m)) = P(m) of the conjecture is valid for all the contents except 
((22), (2),0) and ((2), (22), (2)) marked with *. 
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